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Abstract 

Unlike minimal SU(5), SO(IO) provides a straightforward path towards gauge 
coupling unification by modifying the renormalization group evolution of the gauge 
couplings above some intermediate scale which may also be related to the seesaw 
mechanism for neutrino masses. Unification can be achieved for several different 
choices of the intermediate gauge group below the SO(IO) breaking scale. In this 
work, we consider in detail the possibility that SO(IO) unification may also provide a 
natural dark matter candidate, stability being guaranteed by a leftover Z 2 symmetry. 
We systematically examine the possible intermediate gauge groups which allow a 
non-degenerate, fermionic, Standard Model singlet dark matter candidate while at 
the same time respecting gauge coupling unification. Our analysis is done at the 
two-loop level. Surprisingly, despite the richness of SO(IO), we find that only two 
models survive the analysis of phenomenological constraints, which include suitable 
neutrino masses, proton decay, and reheating. 



1 Introduction 


One of the often quoted motivations for supersymmetry (SUSY) is its ability to improve 
the possibility for gauge coupling unification at the grand unified (GUT) scale, which 
is not possible in minimal SU(5) [1]. However, SO(IO) has the built-in possibility for 
achieving gauge coupling unification through several potential intermediate-scale gauge 
groups [2-4], Of course, low-energy SUSY has many other motivations including the 
presence of a dark matter (DM) candidate [5] whose stability is insured if -R-parity is 
conserved. However, under very generic conditions, non-SUSY SO(IO) models also possess 
a remnant Z 2 symmetry when an intermediate-scale U(l) symmetry is broken [6-11]. 
Thus, several modest extensions of minimal SO(IO) may also allow for the possibility of 
DM. 

In building a successful SO(IO), we must also require that the GUT and intermediate 
mass scales be sufficiently large so as to ensure a proton lifetime and neutrino masses 
compatible with experiment. Unfortunately, these requirements are not realized for ev¬ 
ery choice of intermediate-scale gauge group. The addition of a new SO(IO) multiplet 
containing a DM candidate will, however, affect the running of the gauge couplings and 
can improve the desired unification of the gauge couplings. For this reason, we suppose 
that the DM candidate be charged under the intermediate gauge symmetries. The cosmo¬ 
logical production of DM could occur, for example, out of equilibrium from the thermal 
bath (nonequilibrium thermal DM or NETDM [4]) in a manner reminiscent of freeze-in 
scenarios [12], This mechanism works with a stable particle which has no interaction with 
the SM particles. Thus, we focus on singlet DM candidates. Further, as scalar DM would 
most assuredly couple to the Standard Model (SM) Higgs, we limit our attention here to 
fermionic DM. 

SO(IO) grand unification is, of course, a general moniker for many candidate theories 
of unification, as there are several possible intermediate gauge groups and several possible 
choices for representations R\ of Higgs fields which break SO(IO) to the intermediate 
gauge group, Gj nt , and then again, several possible choices of representations R 2 for the 
Higgs fields which break G int down to the SM. Furthermore, there are several possible 
choices for the representation which contains DM. Thus, it may seem that DM in SO(IO) 
models is a rather robust and generic feature. However, if we insist on maintaining gauge 
coupling unification at a suitably high scale to guarantee proton stability, the number of 
models is dramatically reduced. In fact, by limiting the dimension of the representation 
containing DM to be no larger than a 210 , we find that only two models survive. 

In this paper, we will systematically examine the possibility for fermionic NETDM in 
SO(10) models, though our conclusions are more general than the specific NETDM model. 
We will discuss the various possible intermediate gauge groups and Higgs representations 
which allow for gauge coupling unification, and we will demonstrate the effect of including 
two-loop running of the renormalization group equations (RGEs). The DM representation 
needs to be split so that only fermions with the appropriate gauge quantum numbers 
survive at low energy. This requires fine-tuning similar to the doublet-triplet separation 
problem in GUTs. We also systematically consider viable DM representation and their 
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effect on the running of the gauge couplings. In all but two distinct models, the presence 
of DM spoils the desired unification of the gauge couplings. 

In the following, we begin by discussing the origin of a discrete symmetry in a variety 
of models with different intermediate gauge groups and the possible representations for 
DM and the splitting of the DM multiplet. In Sec. 3, we first demonstrate gauge coupling 
unification in these models (without DM) and show the effect of including the two-loop 
functions in the RGE running and one-loop threshold effects. We next consider the 
question of gauge coupling unification in the presence of a DM multiplet. In Sec. 4, we 
discuss the criteria which select only two possible models in a specific example of the 
NETDM scenario [4]. The phenomenological aspects of these models including neutrino 
masses, proton decay, and the production of DM through reheating after inflation will be 
discussed in Sec. 5. We also consider the case where the DM field is a singlet under the 
intermediate gauge groups in Sec. 6. Our conclusions will be given in Sec. 7. 

2 Candidates 

We assume that the SO(IO) gauge group is spontaneously broken to an intermediate 
subgroup G'int at the GUT scale Mqut, and subsequently broken to the SM gauge group 
Gsm at an intermediate scale M int : 

SO(IO) —> Gi n t —* Gsm <8> ZW , (1) 

with Gsm = SU(3) c ®SU(2) L ®U(l) y . The Higgs multiplets which break SO(IO) and G int 
are called Ri and R 2 , respectively. In addition, we require that there be a remnant discrete 
symmetry Z^v that is capable of rendering a SM singlet field to be stable and hence account 
for the DM in the Universe [10,11]. The mechanism for ensuring a remnant Z^v is discussed 
in detail in Sec. 2.1, and the possible intermediate gauge groups that accommodate the 
condition are summarized in Sec. 2.2. 

If, moreover, the DM couplings are such that the candidate is not in thermal equilib¬ 
rium at early times, as in the NETDM scenario, we obtain stringent constraints on the 
model structure. We will consider this subject in Sec. 2.3. 

2.1 Discrete symmetry in SO(IO) 

SO(IO) is a rank-five group and has an extra U(l) symmetry beyond U(l)y in the SM 
gauge group. The U(l) charge assignment for fields in an SO(IO) multiplet is determined 
uniquely up to an overall factor. We define the normalization factor such that all of the 
fields (f)i in a given model have integer charges Qi with a minimum non-zero value of 
\Qi\ equal to +1. Now, let us suppose that a Higgs field (j) H has a non-zero charge Qh- 
Then, if Qh — 0 (mod. N) with N > 2 an integer, the U(l) symmetry is broken to a 
Z N symmetry after the Higgs field obtains a vacuum expectation value (VEV) [7-9]. One 
can easily show this by noting that both the Lagrangian and the VEV {(J>h) are invariant 


2 



under the following transformations: 


((p H ) ->■ exp ( t 2 n ® H \(j) H ) = {<Ph) ■ 


( 2 ) 


4>i -> exp 


i2nQi \ 


4*i 


Thus, an SO(IO) GUT may account for the stability of DM in terms of the remnant Zjy 
symmetry originating from the extra U(l) gauge symmetry. 

The next task is to determine which type of irreducible representations for the Higgs 
held /fin can be exploited to realize the discrete symmetry. To that end, we follow the 
discussion presented in Ref. [13]. The discussion is based on the Dynkin formalism of the 
Lie algebra [14], 1 Since the rank of SO(IO) is five, we have five independent generators 
which can be diagonalized simultaneously. We denote them by Hi (i = 1,..., 5). They 
form the Cartan subalgebra of SO(IO). Each component of a multiplet is characterized 
by a set of eigenvalues of the generators, /tj (i = 1,..., 5), called weights. We also define 
the weight vector = (/ii,... ,yU 5 ). The weights in the adjoint representation are called 
roots ccj, with a = (cti,... , 0 : 5 ) the root vector. Among the root vectors, a set of Eve 
linearly independent vectors play an important role. They are called simple roots, cx. t 
(?’ = 1,... 5), and expressed by the Dynkin diagrams. In what follows, we consider the 
weight and root vectors in the so-called Dynkin basis. In this particularly useful basis, a 
weight vector /1 is expressed in terms of a set of Dynkin labels given by 


2 a-i ■ n 



( 3 ) 


It turns out that the Dynkin labels are always integers. For example, the highest weight 
of the 16 in SO(IO) is expressed as (0 0 0 0 1), while that of the 10 is given by (1 0 0 0 0). 

On the other hand, it is convenient to express the Cartan generators Hj in the dual 
basis, where they are expressed in terms of five-dimensional vectors [hi 1 ,..., hi 5 ] such that 
their eigenvalues for a state corresponding to the weight /1 are given by 


5 

Hi{p) = hijhj ■ 

3 =1 


( 4 ) 


We choose the five linearly independent Cartan generators as follows: 

#1 = I[i 2 211], 

H>- [(On ill], 

Hi — 1[—2 0 3 - 1 1], 
o 

H 5 = [2 0 2 1 - 1] . (5) 

1 For a review and references, see Refs. [15,16]. We follow the convention of Ref. [15] in this paper. 
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Here, Hi and H 2 correspond to the SU(3)c Cartan generators A 3 /2 and A 8 / 2 , respectively, 
where A a (A = 1,..., 8 ) are the Gell-Mann matrices; 77 3 and H 4 are the weak isospin 
and hypercharge, T 3 l and Y, respectively . 2 H 5 is related to the B — L charge as H 5 = 
—5 (B — L) +4 Y. The additional U(l) symmetry required to generate a discrete symmetry 
is provided by a linear combination of the Cartan generators containing H 5 . Following 
Ref. [13] (see also Ref. [ 8 ]), we define the extra U(l) charge Q 1 by 

Qi = —H 4 -\h 5 = [0 0 - 1 0 0], (6) 

This U(l) charge can be also written as Qi = (B — L) — 2 Y. One can readily find that 
all of the components in 10 and 16 have the U(l) charges of either 0 or ±1. 

Now we consider possible representations, R 2 , for 0# discussed above. First, let us 
determine the possible weight vectors corresponding to the component of 0 # that can 
have a VEV without breaking the SM gauge group. Namely, such a component has a zero 
eigenvalue for Hi (% = 1,...,4). This condition tells us that the corresponding weight 
vectors have the following form: 

\i N = (—TV N — N 0 N) . (7) 


The Q\ charges of the vectors are then given by 

Qi{.Hn) = N . 


( 8 ) 


It is found that the smallest irreducible representation that contains the weight vector 
li n has the highest weight 3 

Ajv = (0 0 0 0 A r ) . (9) 

Its dimension is 16 , 126 , 672 , ... for N = 1 , 2, 3,. .. ■, respectively . 4 To obtain a 
symmetry, IV > 2 is required. Thus, as long as we consider relatively small representations 
(such as those with dimensions not exceeding 210 ), 126 is the only candidate 5 for the 
representation of 0#. In this case, the remnant discrete symmetry is Z 2 . 6 

Under the Z 2 symmetry, the SM left-handed fermions are even, while the SM right- 
handed fermions as well as the Higgs held are odd. One can easily show that this symmetry 
is related to the product of matter parity Pm = (— 1 ) 3 ( B ~U [19] an d the U(l)y rotation by 
67 T, e 6 * 71-5 . Thus, if a SM-singlet fermion (boson) has an even (odd) parity, the remnant Z 2 
symmetry makes the particle stable. In Table 1, we summarize irreducible representations 
that contain fijy. We only show those that have dimensions less than or equal to 210. 
From the table, we find that a singlet fermion in a 45 , 54 , 126 , or 210 representation, or 
a singlet scalar boson in a 16 or 144 representation, can be a DM candidate. 

2 In the case of the flipped SU(5) scenario [17,18], the weak hypercharge is given by Y = —^(H 4 + H 5 ). 

3 In fact, we obtain pjv by subtracting the root vector (1 — 1 1 0 0) from Ajv JV times. 

4 The dimension of An for any N is given by 


dim(AAr) 


(1 + N) 1 + 







( 10 ) 
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Table 1: Irreducible representations containing Miv- 



Representation 

Highest weight 

Z 2 

Mo 

45 

(0 1 0 0 0) 

+ 


54 

(2 0 0 0 0) 

+ 


210 

(0 0 0 1 1) 

+ 

Mi 

16 

(0 0 0 0 1) 

— 


144 

(10 0 10) 

— 

M2 

126 

(0 0 0 0 2) 

+ 


Note that although we need a 126 Higgs held to break the extra U(l) symmetry 
and produce a remnant Z 2 symmetry, other Z, 2 -even singlet fields, 45, 54, 210, etc., can 
have VEVs simultaneously without breaking the Z 2 symmetry. While the latter do not 
break the Z 2 symmetry, as discussed above, they are not capable of producing it, thus 
requiring the 126. We will use such fields to obtain an adequate mass spectrum and a 
non-degenerate DM candidate, as discussed in Sec. 2.3 and Sec. 4. R -2 will therefore refer 
to all representations at the intermediate-scale which are responsible for either symmetry 
breaking or intermediate scale masses and may be a combination of the 126 and other 
representations listed in Table 1 with positive Z 2 charge. 

2.2 Intermediate gauge group 

As shown in Eq. (1), the extra U(l) symmetry is assumed to be broken at the intermediate 
scale, i.e., the 126 Higgs held acquires a VEV of the order of M illt . Thus, the intermediate 
gauge group G- mt should be of rank hve. In Table 2, we summarize the rank-hve subgroups 
of SO(IO) and the Higgs multiplets R\ whose VEVs break SO(IO) into the subgroups. 
Again, we only consider the representations whose dimensions are less than or equal to 
210. Here D denotes the so-called D-parity [20], that is, a Z 2 symmetry with respect to 
the exchange of SU(2) L -B- SU(2) R . D-parity can be related to an element of SO(IO) [20] 
under which a fermion held transforms into its charge conjugate. In cases where the 
D-parity is not broken by R .\, it is subsequently broken by R 2 at the scale of M int . In 
the NETDM scenario, the reheating temperature is always below Mmt, and therefore any 
cosmological relics [6] due to the breaking of D-parity will be harmless. Note that the 
VEVs of the Ri Higgs helds are even under the Z 2 symmetry considered in Sec. 2.1. Thus, 
there is no danger for this Z 2 symmetry to be spontaneously broken by the Ri Higgs helds. 


5 The next-to-smallest representation including p ,2 is 1728 with the highest weight (10 0 1 1). 
6 For earlier work on the remnant Z 2 symmetry in SO(10), see Ref. [6]. 
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Table 2: Candidates for the intermediate gauge group Gi nt . 


Gint R\ 

SU(4) C ®SU(2) L ®SU(2) R 210 

SU(4)c 0 SU(2 ) l ® SU( 2 ) r 0 D 54 

SU(4) C 0 SU(2) l 0 U(l) fl 45 

SU(3) c ®SU(2) L ®SU(2) R ®U(l) s _ L 45 

SU(3) C 0 SU(2) L 0 SU(2) R 0 U(1 ) B -l 0 D 210 
SU(3) C 0 SU(2) L 0 U(l) fl 0 U(1)b-l 45, 210 

SU(5) 0 U(l) 45, 210 

Flipped SU(5) 0 U(l) 45, 210 


2.3 Fermion dark matter and degeneracy problem 

In the NETDM scenario, the DM should not be in thermal equilibrium. This requirement 
disfavors scalar DM candidates, since a scalar, <f >, can always have a quartic coupling with 
the SM Higgs held H: ^<ph\4>\ 2 \ H\ 2 . Unless |A^| is extremely small for some reason, 
this coupling keeps scalar DM in thermal equilibrium even when the temperature of the 
Universe becomes much lower than the reheating temperature. Therefore, we focus on 
fermionic DM in this paper. Following the discussion in Sec. 2.1, the DM candidate should 
be contained in either a 45, 54, 126, or 210 representation. 

Below the GUT scale, components in an SO(IO) multiplet can obtain different masses. 
We assume that only a part of an SO(IO) multiplet which contains the DM candidate 
and forms a representation under G illt has a mass much lighter than the GUT scale. We 
denote this representation by Rdm- Such a mass splitting can be realized by the Yukawa 
coupling of the DM multiplet with the Ri Higgs held. After the Ri Higgs obtains a VEV, 
the Yukawa coupling leads to an additional mass term for the SO(IO) multiplet, which 
gives different masses among the components. By carefully choosing the parameters in 
the Lagrangian, we can make only Rdm light. This will be discussed in detail in Sec. 4. 

As will be seen in Sec. 3.1, without Rdm, SO(IO) GUTs often predict a low value of 
either Mqdt or M int , which could be problematic for proton decay or the explanation 
of light neutrino masses, respectively. In order to affect the RGE running of the gauge 
couplings and possibly increase the mass scales for both M int and Mqdt, the DM should 
be charged under Gj nt . In Table 3, we summarize possible candidates for Rdm for each 
intermediate gauge group. Above the intermediate scale, all of the components have an 
identical mass. In fact, it turns out that the degeneracy is not resolved at tree level even 
after the intermediate gauge symmetry is broken. This is because the SO(IO) multiplets 
which contain Rdm displayed in the table cannot have Yukawa couplings with the 126 
Higgs; such a coupling is forbidden by the SO(IO) symmetry. Thus, the effects of sym- 
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Table 3: Candidates for the NETDM. 


Dint. 

7?dm 

SO(10) 

SU(4)c ® SU(2 )l ® SU(2 ) R 

( 1 , 1 , 3 ) 

45 


( 15 , 1 , 1 ) 

45 , 210 


( 10 , 1 , 3 ) 

126 


( 15 , 1 , 3 ) 

210 

SU(4) C ® SU(2)i <8 ) U(1)b 

( 15 , 1 , 0 ) 

45 , 210 


( 10 , 1 , 1 ) 

126 

SU(3) C <g> SU(2) L ® SU(2 ) r ® U(1) B _ L 

( 1 , 1 , 3 , 0 ) 

45 , 210 


( 1 , 1 , 3 , - 2 ) 

126 

SU(3) C <g> SU(2) L ® U(l) fl <g> U(1 ) B -l 

( 1 , 1 , 1 , -2) 

126 

SU(5) ® U(l) 

( 24 , 0 ) 

45 , 54 , 210 


( 1 , -10) 

126 


( 75 , 0 ) 

210 

Flipped SU(5) ® U(l) 

( 24 , 0 ) 

45 , 54 , 210 


( 50 , -2) 

126 


( 75 , 0 ) 

210 


metry breaking by the 126 Higgs VEV cannot be transmitted to the mass of the Rbm 
mnltiplet at tree level, and a simple realization of DM in Rbm makes its components 
degenerate in mass. 

Such a degenerate mass spectrum is problematic. Since the degenerate mnltiplet 
contains particles charged under the SU(3)c<8>U(1)em gauge group, they will be in thermal 
equilibrium. In general, these components have quite a long lifetime, and thus their 
thermal relic density conflicts with various observations. To see this, let us consider the 
( 1 , 1 , 3 ) Dirac fermion multiplet (^°, V ,=b ) hi the SU(4)o <g) SU(2)z, <S) SU(2 )theory, which 
originates from the 45 representation of SO(IO), as an example. As mentioned above, 
they have an identical mass M at tree level, and the mass difference AM induced by the 
radiative corrections can be estimated as 

AM~-Mln( —) - O.OlxM, (11) 

4vr \ M J ’ y J 

where an is the U(l) gauge fine-structure constant. The charged components '0 ± can 
decay into the neutral DM if 0 only through the exchange of the intermediate-scale gauge 
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f 


Figure 1: Diagram responsible for the decay of - 0 " 
SM particles. 


into the DM if 0 , f and f denote the 


bosons as shown in Fig. 1. We estimate the decay width as 


r(i/> + </>"//') ~ 


a 


R 


(A Mf 


n 


M w K 


( 12 ) 


where an = and gn and M Wr are the coupling and the mass of the intermediate 

gauge boson Wr, respectively. Then, for example, when the DM mass is 0(1) TeV and 
the intermediate scale is 0(1O 13 ) GeV, the lifetime of is much longer than the age of 
the Universe, and thus cosmologically stable. The abundance of such a stable charged 
particle is stringently constrained by the null results of the search for heavy hydrogen 
in sea water [21], The DM multiplcts in other cases may also be accompanied by stable 
colored particles, whose abundance is severely restricted as well. If the intermediate 
scale is relatively low, the charged/colored particle can have a shorter lifetime. Even in 
this case, their thermal relic abundance should be extremely small in order not to spoil 
the success of the Big-Bang Nucleosynthesis (BBN). Quite generally, a degenerate mass 
spectrum leads to disastrous consequences. We refer to this problem as the “degeneracy 
problem’’ in what follows. 

To avoid the degeneracy problem, we need to make the charged/colored components 
heavy enough so that they are not in thermal equilibrium and have very short lifetimes. 
To that end, it is natural to explore a way to give them masses of 0(M int ) by using the 
effects of the intermediate symmetry breaking. There are several solutions. One of the 
simplest ways is to introduce an additional Higgs held that has a VEV of the order of 
M int . For this purpose, we can use a 45 , 54 , or 210 held, as discussed in Sec. 2.1. The 
Yukawa coupling between the Higgs and the DM then yields the desired mass splitting. By 
hne-tuning the coupling we can force only the DM to have a mass much below M- mt while 
the other components remain around the intermediate scale'. Though other mechanisms 
are possible, we adopt this approach in this work. Concrete realizations of the mechanism 
are illustrated in Sec. 4. 

7 This fine-tuning is similar (though somewhat less severe) to the fine-tuning associated with the 
doublet-triplet separation to insure a weak scale Higgs boson. 
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Another solution to the degeneracy problem involves the use of higher-dimensional 
operators that include at least two 126 fields. One would expect that such operators 
suppressed by the Planck scale, Mpi, always exist. These Planck-suppressed operators 
can give rise to a mass difference of 0(M^ nt /A / Ipi). Another mechanism to generate higher¬ 
dimensional operators is to introduce a vector-like fermion which has a Yukawa coupling 
with the DM and the 126 Higgs. By integrating out the fermion, we obtain dimension- 
five operators which give an 0(M^ nt /M{ ei ) mass difference, where Mf er is the mass of the 
additional fermion. Moreover, the higher-dimensional operators can be induced at the 
loop level, which gives rise to a (D(aG\jTM? nt / (47tMgut)) mass difference, where «gut = 
5 , Gu r r/(4 7r ) is the fine-structure constant of the unified gauge coupling g G UT- Realization 
of these scenarios will be discussed elsewhere. 


3 Gauge coupling unification 

As is well known, gauge coupling unification can be realized in SO(IO) GUTs with an 
intermediate scale [ 2], 8 Once the intermediate gauge group as well as the low-energy 
matter content is given, one can determine both the intermediate and GUT scales by 
requiring gauge coupling unification. In what follows, we reevaluate these scales in the 
SO(IO) GUT scenarios with different intermediate gauge groups and up-to-date values for 
the input parameters. Then, in Sec. 3.2, we study the effects of the DM and the interme¬ 
diate Higgs multiplets on gauge coupling unification. We will find that the requirement 
of gauge coupling unification severely constrains the NETDM models. 

3.1 Gauge coupling unification with the intermediate scale 

To begin with, let us briefly review SO(IO) GUTs with an intermediate gauge group. In 
SO(IO) GUTs, the SM fermions as well as three right-handed neutrinos are embedded into 
three copies of the 16 spinor representations, while the SM Higgs boson is usually included 
in a 10 representation. At the GUT scale, the SO(IO) GUT group is spontaneously broken 
into an intermediate gauge group. Subsequently, the intermediate Higgs multiplet breaks 
it into the SM gauge group at the intermediate scale. In the following analysis, we work 
with the so-called extended survival hypothesis [23,24]; that is, we assume that a minimal 
set of Higgs multiplets necessary to realize the symmetry breaking exists in low-energy 
region. Above the intermediate scale, the presence of the additional Higgs multiplet and 
intermediate gauge bosons change the gauge coupling running from that in the SM. This 
makes it possible to realize gauge coupling unification in this scenario. 

As displayed in Table 2 , the intermediate gauge groups relevant to our discussion are 
divided into two classes; those which contain the SU(5) group as a subgroup, and those 
which do not. The former class is, however, found to be less promising. In the case of 
ordinary SU(5) (8) U(l), the SM gauge couplings should meet at the intermediate scale, 
though they do not, as is well known. Failure of gauge coupling unification is also found in 

8 For a review, see Refs. [3,22]. 
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the flipped SU(5) case. This conclusion cannot be changed even if one adds the DM and 
Higgs mnltiplets in the case of ordinary SU(5). In the flipped SU(5) case, the addition of 
the DM and Higgs mnltiplets may yield gauge coupling unification. However, it turns out 
that the intermediate mass scale is as high as (9(10 1 ') GeV in such cases. Since the masses 
of the right-handed neutrinos are expected to be if M int = 0(H) 1 '), the simple 

seesaw mechanism [25] cannot explain the neutrino masses required from the observation 
of the neutrino oscillations. However, the GUT scale tends to be close to the Planck scale, 
and one may need to rely on a double seesaw to explain neutrino masses [18,26]. We do 
not consider these possibilities in the following discussion. 

The other class of the intermediate gauge groups is related to the Pati-Salam gauge 
group [27]. Therefore, it is useful to decompose the SO(IO) multiplets into mnltiplets of 
the SU(4)c <E> SU(2)l 8 SU(2)r gauge group. The 16 spinor representation in SO(IO) is 
decomposed into a (4, 2,1) and (4,1, 2) of SU(4)c 8 SU(2)^ (8) SU(2 )r. We denote them 
by and T'p, respectively, in which the SM fermions are embedded as follows: 


4 4 4 
4 4 4 e L J ’ 




f 4i 

\~ U R1 


( 1 C R2 

a, 

CO 

4 

~ u< k-2 

U R3 

-4 


(13) 


where the indices represent the SU(3)c color and C indicates charge conjugation. The 
SM Higgs field is, on the other hand, embedded in the (1,2,2) component of the ten¬ 
dimensional representation. As discussed in Ref. [28], to obtain the viable Yukawa sector, 9 
we need to consider a complex scalar 10c for the representation, not a real one. Thus, 
(1,2,2) is also a complex scalar multiplct and includes the two Higgs doublets. In the 
following calculation, we regard one of these doublets as the SM Higgs boson, and the other 
is assumed to have a mass around the intermediate scale. The SU(4)c® SU(2)c® SU(2) R 
gauge group is broken by the VEV of the (10,1,3) component in the 126c- In the 
presence of the left-right symmetry, we also have a (10, 3,1) above the intermediate scale. 
We assume that the (10, 3,1) field does not acquire a VEV, with which the constraint 
coming from the p-parameter is avoided. From these charge assignments, one can readily 
obtain the quantum numbers for the corresponding fields in the other intermediate gauge 
groups, since they are subgroups of the SU(4)c 8 SU(2)c 8) SU(2)#. 

With this field content, we study whether the gauge coupling unification is actually 
achieved or not for the first six intermediate gauge groups listed in Table 2. We perform 
the analysis by using the two-loop RGEs, which are given in Appendix B. We will work 
in the DR scheme [30], as there is no constant term in the intermediate and GUT scale 
matching conditions. The input parameters we use in our analysis are listed in Table 7 in 
Appendix A. By solving the RGEs and assuming gauge coupling unification, we determine 
the intermediate scale M int , the GUT scale Mqut, and the unified gauge coupling constant 
fl'GUT- If we fail to find the appropriate values for these quantities, we will conclude that 
gauge coupling unification is not realized in this case. To determine their central values 

9 For a general discussion on the Yukawa sector in SO(IO) GUTs, see Refs. [28,29]. 
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Tabic 4: log 10 (M int ), log 10 (Mgut), and gcur- For each Gi nt , the upper shaded (lower) 
row shows the 2-loop (1-loop) result. M int and Mqut are given in GeV. The blank entries 
indicate that gauge coupling unification is not achieved. 


Gint 


logio(M in t) 

logic (Mqut) 

<?GUT 

SU(4) c g)SU(2) L 0SU(2) H 


11.17(1) 
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0.52738(4) 
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SU(4) C g SU(2) l g SU(2) fi 

0h 
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14.95(1) 
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15.23(3) 

0.5520(1) 

SU(4) c 0SU(2) l 0U(1) k 


11.35(2) 

14.42(1) 

0.5359(1) 



11.23(1) 

14.638(8) 

0.53227(7) 

SU(3) c 0SU(2) l 0SU(2) h 

® U(1 )b-l 

9.46(2) 

16.20(2) 

0.52612(8) 



8.993(3) 

16.68(4) 

0.52124(3) 

SU(3) c 0SU(2) l 0SU(2) h 

0 U(1)b-l <0 D 

10.51(1) 

15.38(2) 

0.53880(3) 



10.090(9) 

15.77(1) 

0.53478(6) 

SU(3) C 0 SU(2 ) l 0 U(l) fl 0 U(l)s_ L 


as well as the error coming from the input parameters, we form a y 2 statistic as 


X 2 = S 


< 2=1 


Wa ~ ^exp ) 2 

a 2 {g 2 ' 


(14) 


a,exp/ 


where g a are the gauge couplings at the electroweak scale obtained by solving the RGEs 
on the above assumption, and g a . eX p are the experimental values of the corresponding 
gauge couplings, with cr(g 2 exp ) denoting their error. The central values of M int , M G ut , 
and gam are corresponding to a point at which y 2 is minimized. 10 

By using the method discussed above, we carry out the analysis and summarize the 
results in Table 4. Here, we show log 10 (M int ), log 10 (McuT), and (?gut- For each inter¬ 
mediate gauge group, the upper shaded (lower) row shows the 2-loop (1-loop) result. 
M int and M GU t are given in GeV. The blank entries indicate that gauge coupling unifi¬ 
cation is not achieved. The uncertainties resulting from the input error are also shown 
in the parentheses. To illustrate our procedure more clearly, we show y 2 as functions of 
log 10 (M int ) (top), log(McuT) (middle), and gcxjT (bottom) in Fig. 2 for two examples of 
intermediate gauge groups. The left panels are for G int = SU(4)<y 0 SU(2)i 0 SU(2)b, 
while the right ones are for G int = SU(4) C 0 SU(2) L g) SU(2) R g D. The y 2 functions for 

10 We also use the y 2 statistics to determine the value of the input Yukawa coupling in a similar manner, 
though it scarcely affects the error estimation of Mqut, and ^gut- 
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Figure 2: \ 2 as functions of \og 10 (M int ) (top), log 10 (M GUT ) (middle), and g gut (bottom). 
Left and right panels are for Gi nt = SU(4)c <S> SJJ(2)l ® SJJ(2)r and Gi nt = SU(4)c <S> 
SU{2)l ® SU{2)r (g> D, respectively. M int and Mgut are given in GeV. 
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Figure 3: Contour plots for the allowed region in the g GUT-\og 10 (M int ), g GUT-\og 10 (M gut) , 
and \og 10 (M G uT)-log 10 (M in t) parameter planes in the top, middle, and bottom panels, 
respectively. Left panels are for Gi nt = SU(A)c ® SU(2 )l <8> SU(2)r, while right ones are 
for G int = SU(A)c ® SU(2 ) l <8) SU{2) R ® D. Stars represent the best-fit point. The colored 
regions correspond to 68, 95, and 99 % C.L. limits determined from A\ 2 — 2.30,5.99,9.21. 
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the other choices of intermediate scale gauge groups will be qualitatively similar. Here 
again, M- mt and Mgut are given in GeV. In each plot, the other two free parameters are 
fixed to their best-fit values. We also plot the one-loop results (shown as dotted curves) 
to show the significance of the two-loop effects. In Fig. 3, we show the y 2 functions pro¬ 
jected down onto 2D planes corresponding to <7GUT-log 10 (M int ), 5 , GUT-log 10 (d^GUT), and 
log 10 (-Mqut) _ log 10 (fHint) in the top, middle, and bottom panels, respectively. Again, the 
left panels are for G'j nt = SU(4)c 8 SU(2)u 8 SU(2)i{, while the ones on the right are for 
Gi n t = SU(4 )c8SU(2)u8SU(2)k 8-D. The stars represent the best-fit point. The uncer¬ 
tainty ellipses represent 68, 95, and 99 % C.L. uncertainties corresponding to Ay 2 = 2.30, 
5.99, and 9.21, respectively. Threshold corrections at M int and Mgut [31] due to the 
non-degeneracy of the particles that have masses of the order of these scales contribute to 
the uncertainties. 11 For a recent discussion of threshold corrections, see Ref. [32], In addi¬ 
tion, we neglect the contribution of Yukawa couplings above the intermediate scale, which 
causes additional error. These are expected to give 0(1)% uncertainty to the results. 

From Table 4, it is found that gauge coupling unification is not achieved in the case of 
G int = SU(3) C 8 SU(2)x, (g) U(1 )r 8) U(1)b-l- Moreover, we find that relatively low GUT 
scales are predicted for G; nt = SU(4)c8SU(2)u8SU(2)^8-D and SU(4)c8SU(2)L8U(l)i?, 
and thus the proton decay constraints may be severe in these cases, as discussed in Sec. 5.2. 
Furthermore, except for Gj nt = SU(4 )q8SU(2)u 8SU(2)^8-D, we obtain low intermediate 
scales, with which it may be difficult to account for the neutrino masses, as explained in 
Sec. 5.1. As we will see below, this situation can be improved in the NETDM models. 

3.2 NETDM and gauge coupling unification 

Next, we look for the NETDM models in which gauge coupling unification is realized with 
an appropriate intermediate unification scale. Here, we require 10 15 < M G ut ;$ 10 18 GeV; 
if Mqut < 10 15 GeV, then proton decays are too rapid to be consistent with proton decay 
experiments, while if Mgut > 10 18 GeV, then gravitational effects cannot be neglected 
anymore and a calculation based on quantum held theories may be invalid around the 
GUT scale. To search for promising candidates, we assume the following conditions. 
Firstly, a model should contain a NETDM candidate shown in Table 3, where only a 
singlet component has a mass much below the intermediate scale. This component does 
not affect the running of the gauge couplings. Secondly, the rest of the components in Rdm 
are assumed to be around M; nt clue to the mass splitting mechanism with an additional 
Higgs multiplet, discussed in Sec. 2.3. At this point, we only assume that there exists an 
extra Higgs multiplet from either the 45, 54 or 210 whose mass is around the intermediate 
scale. Whether the VEV of the extra Higgs actually gives rise to the mass splitting or 
not will be discussed in the subsequent section. Thirdly, we require that only the SM 
fields, the intermediate gauge bosons, Rdm, and R 2 are present below the GUT scale. For 
example, if we consider the (1, 1, 3) DM of the 45 given in the Erst column in Table 3, 
then we suppose that all of the components of the 45 except Rdm = (1,1,3) should have 

11 Note that the intermediate scale in the left-right symmetric theories does not depend on physics 
beyond as discussed in Appendix C. 
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Tabic 5: Models that realize the gauge coupling unification. M int and Mqut are given in 
GeV. All of the values listed here are evaluated at one-loop level. 


SU (4) g G>SU (2) L<g)SU (2)r 

Rdm 

R 2 

logic (Mint) 

l°gl 0 (M G UT) 

9GVT 

(1,1, 3 ) w 

(10,1, 3) g 
( 1,1, 3)j* 

10.8 

15.9 

0.53 

(1,1,3)23 

(10,1, 3) g 
( l,l,3) fl 

9.8 

15.7 

0.53 

SU(4) c ®SU(2) L ®SU(2) fl <g> D 

Rdm 

r 2 

logic (Mint) 

l°glo(M GU T) 

S'GUT 


(10.1.3) C 




(15,1,1 ) w 

(10, 3, l)c 
(15,1,1)22 

13.7 

16.2 

0.56 


(10.1.3),; 

(10, 3,1) G 




(15,l,l) w 

(15,1,3)22 

(15,3,1)22 

14.2 

15.5 

0.56 


(10,l,3) c 

(10,3,1V 




(15,1, 1)d 

(15,1,3)22 

(15,3,1)22 

14.4 

16.3 

0.58 

SU(3)c®SU(2) l ®SU(2) r ®U(1) b l 

-Rdm 

r 2 

logio(M int ) 

1 0 Sio(^gut) 

5Gut 

(1,1, 3, 0 )w 

(1,1,3,-2V 

(1,1,3,0)22 

6.1 

16.6 

0.52 


masses around the GUT scale. This condition is corresponding to the requirement of the 
minimal fine-tunings in the scalar potential to realize an adequate mass spectrum. 

With these conditions, we then search for possible candidates by using the one-loop 
analytic formula given in Appendix C. In Table 5, we summarize the held contents that 
satisfy the above requirements, as well as the values of log 10 (M int ), log 10 (M GU T), and 
S'gut, with M int and M GU t hi GeV. All of the values are evaluated at one-loop level. 
Here the subscript II. C. W, or 1) of each multiplet indicates that it is a real scalar, a 
complex scalar, a Weyl fermion, or a Dirac fermion, respectively. As for the intermediate 
Higgs fields, R 2 , listed in Table 5, (10,1,3V and (1,1,3, —2) G are from the 126 Higgs 
held, while all other representations included in R 2 are extra Higgs helds introduced to 
resolve the degeneracy problem. For the additional Higgs helds, we only show the real 
scalar cases for brevity. Indeed, we can also consider complex scalars for the Higgs helds 
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and find that gauge coupling unification is also realized in these cases, where both the 
intermediate and GUT scales are only slightly modified. 


4 Models 

In the previous section, we have reduced the possible candidates to those presented in 
Table 5. In this section, we study if any of those models are viable; i.e., we check if they 
actually offer an appropriate mass spectrum to realize the NETDM scenario, with the 
charged/colored components in R G m acquiring masses of O (M int ). 

First, let us consider the ( 1 , 1,3 ) w/d DM representation in the SU(4)c <E> SU(2) L <g) 
SU(2)/? gauge theory. To split the masses in the (1,1, 3) multiplet ^ r , we need to couple 
the DM with the (1, 1,3)/? Higgs 0 r , with r denoting the SU(2)/? index. Since the fields 
transform as triplets under the SU(2)/? transformations, to construct an invariant term 
from the fields, the indices should be contracted anti-symmetrically; he., the coupling 
should have a form like 

e pqr (YYY q (Y . (15) 

Then, if Y r is a Majorana fermion, the above term always vanishes. Thus, Y r should be 
a Dirac fermion, that is, (1, 1,3)// is the unique candidate for NETDM in this case. 

Next, we study the terms in the SO(IO) Lagrangian relevant to the masses of the fields 
much lighter than the GUT scale. In SO(IO), (1,1, 3)//, (1,1, 3 )/?, and (10,1, 3)c are in¬ 
cluded in the 45//, 45 /?, and 126c, respectively. The SO(IO) gauge group is spontaneously 
broken by the 210/? Higgs held ( Ri ) into the SU(4) C ® SU(2)/, <g) SU(2)/? intermediate 
gauge group. As is usually done in the intermediate scale scenario, we fine-tune the Higgs 
potential so that the (1,1,3)/? and (10, l,3)c Higgs fields have masses around the in¬ 
termediate scale. This can always be performed by using the couplings of the 45/? and 
126c fields with the 210/?. Higgs field, which acquires a VEV of the order of the GUT 
scale. Similarly, we give desirable masses to the fields in (1,1,3)// by carefully choosing 
the couplings of the 45// fermion with the 45/? and 126c Higgs fields. Here, the relevant 
interactions are 

U int = —M 45 d 45^45// - fy 45 45//45//45/? - y 210 45^45//21O/? . (16) 

Notice that 45// does not couple to the 126c field, as already mentioned in Sec. 2.3. After 
the Ri = 210/? Higgs field gets a VEV (210/?) = u 2 io, the interactions in Eq. (16) lead to 
the following terms: 12 


Ant , (17) 

with Mdm = M 45d + 1 / 210 ^ 210 /^. Here, ?// and 0 r denote the (1,1,3)// and (1,1,3)/? 
components in 45// and 45/?, respectively. We find that although M 45o and n 2 io are 
expected to be 0(M G ut), we can let M dm be much lighter than the GUT scale by carefully 
choosing the above parameters so that they cancel each other. In addition, it turns out 

12 For the computation of the Clebsch-Gordan coefficients, we have used the results given in Ref. [33]. 


16 



that the mass term of the (1,3, 1)# component in 45# is given by M 45d — y 2 io^io/vA 
Thus, even if we fine-tune M 45d and 1/210 to realize Mdm Msut, the mass of (1,3,1)23 
is still around the GUT scale. This observation reflects the violation of the D-parity 
in this model. At this point, all of the components in ip r have identical masses (the 
“degeneracy problem”). Once the neutral component of 0 r acquires a VEV (0 3 ) = u 45 , 
which is assumed to be 0(M int ), the second term in Eq. (17) gives rise to additional mass 
terms for t/k. These are 


Ant -MdmV’V 0 - M + ^ + - M^MA , (18) 


where M± = M dm =|= y 45 u 45 , and and ■/A are the neutral and charged components, 
respectively. 13 The above expression shows that the VEV of the 45# Higgs held indeed 
solves the degeneracy problem; if Mdm M int and y 4 5U45 = 0(M int ), then the charged 
components acquire masses of while the neutral component has a mass much 

lighter than M int . Thus, we obtain the mass spectrum we have assumed in the previous 
section. 

In the next example, we consider the DM representation Rbm = (15,1,1 )w with R 2 = 
(10,1, 3)c©(10, 3, l)c©(15,1, 1)# in the left-right symmetric SU(4)c®SU(2)#0SU(2)# 
gauge theory. In this case, R± = 54#. We assume that the (15,1,1 )w is a part of the 
45w, while both (10, l,3)c and (10,3,1 )c are part of the 126 c . The couplings of the 
DM with the Higgs fields, as well as its mass term, are then given by 

Ant = - 45 ir 45 ir - ^ 45 w 45 w 54 r - ^45^45^210# + h.c. , (19) 


Here, (15,1,1)/? is included in the 210# held; we cannot use a 45# in this case since 
the Weyl fermion 45 w has no coupling to the 45#. 14 As before, below the GUT scale, 
the VEV of 54#, u 54 , gives a common mass M to the (15,1,1 )w multiplet with M = 
M^ w ~~ 2 / 54 ^ 54 /\/l5- We can take M = by fine-tuning M 45vk and y 54 U 54 . The 

above Lagrangian then reduces to 

Ant ->■ ~y^ A ^ A + + h.c. , (20) 


where i\) A and (p A denote the (15,l,l)w and (15,1,1)# helds, respectively, with if} = 
i\) a T a and (f) = <p A T A ; A,B,C = 1,... 15 are the SU(4) adjoint indices and T A are the 
SU(4) generators. The mass degeneracy in this case is resolved by the VEV of the 210# 
held, 


(A 


^diag(l,l,l,-3), 


( 21 ) 


with which Eq. (20) leads to 


£in. -> - -V’V - ^ 9 V ~ + h.c. , (22) 

13 Note that since ip r are Dirac fermions, (t/d) C 7 ^ V ’ 0 an< i (tA) C A A 

14 It is also possible to embed (15,1, l) w into 210 w and (15,1,1)# into 45#. The phenomenology in 
this case is the same as that discussed in the text. 
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where 0 °, g A , £“, and £ a are the color singlet, octet, triplet, and anti-triplet components 
in (15,1,1 )w, respectively, with a denoting the color index, and 


y /2 

Mdm — M H— —IJ210V210 , 

(23) 

Mg — M 3 ^ 2 / 210^210 , 

(24) 

1 

Me — M + ^U 2 io v 2 w ■ 

(25) 

Therefore, by carefully adjusting 2 / 210 ^ 210 , we can make the DM -0° 

much lighter than M int 


while keeping the other components around the intermediate scale. 

There are two more possible representations for Rbm for the left-right symmetric 
SU(4)c<g)SU(2)#Cg>SU(2)# intermediate gauge group given in Table 5, namely (15,1,1 ) w / d - 
In this case, however, one can readily conclude that the degeneracy problem cannot be 
solved by the (15,1,3)# and (15,3,1)# Higgs fields. This is because the Yukawa cou¬ 
plings between the DM and these Higgs fields are forbidden by the intermediate gauge 
symmetry. As a consequence, we can safely neglect these possibilities. 

Finally, we discuss the model presented in the last column in Table 5. We again find 
that the (1,1,3,0)# Higgs held does not yield a mass difference among the components 
in the (1,1,3, 0)^ DM multiplet, since the operator in Eq. (15) vanishes when the DM 
is a Weyl fermion. Thus, we do not consider this model in the following discussion. 

As a result, we obtain two distinct models for NETDM within SO(IO). We summarize 
these two models in Table 6 . We call them Model I and II in what follows. Here, M mt and 
Mgut are given in GeV, and all of the values are evaluated with two-loop RGEs and differ 
somewhat from the one-loop values given in Table 5. The errors shown in the parentheses 
arise from uncertainties in the input parameters. In addition, we again expect threshold 
corrections at M int and Mgut- Furthermore, we neglect the contribution of Yukawa 
couplings to the RGEs above the intermediate scale, and this also will contribute to the 
theoretical error. We estimate that these two sources cause 0(1)% uncertainties in the 
values displayed in Table 6 . From these results, we find that the presence of the DM 
component as well as the extra Higgs bosons can significantly alter the intermediate and 
GUT scales , 15 because of their effects on the gauge coupling running. To illustrate this 
more clearly, in Fig. 4 we show the running of gauge couplings in each theory. The left 
and right panels of Fig. 4 correspond to Model I and II, respectively. In each figure, solid 
(dashed) lines show the case with (without) DM and additional Higgs bosons. The blue, 
green, and red lines represent the running of the U(l), SU(2) and SU(3) gauge couplings, 
respectively, where the U(l) fine-structure constant aq is defined by 


1 


3 1 


2 1 


5 CK 2 # 5 «4 ’ 


(26) 


15 However, their existence hardly changes the intermediate scale in Model II, which is clarified in 
Appendix C. 
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while the SU(3)c coupling a 3 is defined by a 3 = a 4 above the intermediate scale. These 
figures clearly show the effects of the extra particles on the gauge coupling running. In 
particular, the GUT scale in Model II is now well above 10 15 GeV, which allows this 
model to evade the proton decay constraints, as will be seen in the subsequent section. 


Table 6: NETDM models. Mi nt and Mgut are given in GeV. All of the values are 
evaluated with the two-loop RGEs. 



Model I 

Model II 

G mt 

SU(4) C <g) SU(2 ) l <8 ) SU(2)^ 

SU(4) C <8 SU(2) l ® SU(2) r ® D 

-Rdm 

(1,1, 3) d m45 D 

(15,1, l)w in 45 w 

Ri 

210 r 

54 r 

R'2 

(10,1, 3)c © (1,1, 3)r 

(10,l,3) c ® (10,3, l) c ® (15,1, 1) R 

loglo(Mnt) 

8.08(1) 

13.664(5) 

logic (M GU t) 

15.645(7) 

15.87(2) 

^GUT 

0.53055(3) 

0.5675(2) 


G irrt = SU(4) c ®SU(2) l ®SU(2) r 



Gint = SU(4)c®SU(2) l ®SU( 2) r ®D 


16 


(a) Model I 



Figure 4: Running of gauge couplings. Solid (dashed) lines show the case with (without) 
DM and additional Higgs bosons. Blue, green, and red lines represent the running of the 
U(l), SU(2) and SU(3) gauge couplings, respectively. 
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5 Phenomenological aspects 

Now that we have obtained the NETDM models, we can study their phenomenological 
aspects and possible implications in future experiments. In Sec. 5.1, we first consider 
whether these models can give appropriate masses for light neutrinos. Next, in Sec. 5.2, 
we evaluate proton lifetimes in each model and discuss the testability in future proton 
decay experiments. Finally, we compute the abundance of DM produced by the NETDM 
mechanism in Sec. 5.3, and predict the reheating temperature after inflation. 


5.1 Neutrino mass 


In SO(IO) GUTs, the Majorana mass terms of the right-handed neutrinos are induced 
after the B — L symmetry is broken. These mass terms are generated from the Yukawa 
couplings of the 16 spinors with the 126 c Higgs held. If the Yukawa couplings are 0( 1), 
then the Majorana mass terms are 0(M int ). On the other hand, in these models, the 
Dirac masses of neutrinos are equal to the up-type quark masses, m u , at the unification 
scale. Therefore, via the seesaw mechanism [25], light neutrino masses are given by 


m: 


m v ~ 




int 


(27) 


In Model II, M int = (9(10 13 ) GeV indeed gives proper values for neutrino masses . 16 How¬ 
ever, in Model I, a low intermediate scale of 0(1O 8 ) GeV yields neutrino masses which 
are too heavy using the standard seesaw expression (27). Thus, Model I is disfavored on 
the basis of small neutrino masses. 

The defect in Model I may be evaded if the (15, 2, 2) component in 126c has a 
sizable mixing with the (1,2,2) Higgs boson and acquires a VEV of the order of the 
electroweak scale. In this case, the neutrino Yukawa couplings can differ from those of the 
up quark, and thus the relation (27) does not hold any more. For sizable mixing to occur, 
the (15, 2, 2) held should he around the intermediate scale. One might think that the 
presence of additional fields below the GUT scale would modify the running of the gauge 
couplings and spoil the above discussion based on gauge coupling unification. However, it 
turns out that both the intermediate and GUT scales are hardly affected by the existence 
of this held, though the unified gauge coupling constant becomes slightly larger. This 
is because its contribution to the one-loop beta function coefficients is Ab 4 = 16/3 and 
Ab 2 L = Ab- 2 R = 5, and thus their difference is very tiny (see the discussion given in 
Appendix C). Therefore, we can take the (15,2,2) to be at the intermediate scale with 
little change in the values of M int and M G ut- The presence of the (15, 2, 2) is also desirable 

16 Note that in a left-right symmetric model such as Model II there is in general also a type-II seesaw 
contribution to m„ from the VEV of the SU(2ffi triplet in the 126(7 . However, we know from constraints 
on the p-parameter that the VEV must be quite small and definitely much smaller than the VEV of 
the SU(2)fl triplet. For example, if the mixing between the SU(2)i and SU(2)^ triplets with the Higgs 
doublets is small, it is safe to assume that the SU(2 )l triplet VEV is small and thus the type-II seesaw 
contribution is subdominant [34]. 
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to account for the down-type quark and charged lepton Yukawa couplings [28,35-37]. In 
addition, the higher-dimensional operators induced above the GUT scale may also affect 
the Yukawa couplings. Constructing a realistic Yukawa sector in these models is saved 
for future work. 

5.2 Proton decay 

Proton decay is a smoking-gun signature of GUTs, and thus a powerful tool for testing 
them. In non-SUSY GUTs, p —y e + n° is the dominant decay mode, which is caused by 
the exchange of GUT-scale gauge bosons. This could be compared with the case of the 
SUSY GUTs; in SUSY GUTs, the color-triplet Higgs exchange usually yields the dominant 
contribution to proton decay, which gives rise to the p —> K + v decay mode [38]. 17 

Since the p —> e + 7r° decay mode is induced by gauge interactions, we can make a robust 
prediction for the partial decay lifetime of this mode. Details of the calculation are given 
in Appendix D. By using the results given there, we evaluate the partial decay lifetime of 
the p — > e + 7r° mode in each theory, and plot it as a function of M X /M G ut ( M x denotes the 
mass of the GUT-scale gauge boson) in Fig. 5. Here, the blue and red solid lines represent 
Models I and II, while the blue and red dashed lines represent the models without the DM 
and extra Higgs multiplets as given in Table 4, namely G'j nt = SU(4)c 0 SU( 2 )£ (8) SU( 2 )j? 
and Gi n t = SU(4 )cG>SU(2)l(8)SU(2)rG)D, respectively. The shaded area shows the region 
which is excluded by the current experimental bound, r(p —> e + 7r°) > 1.4 x 10 34 years 
[40,41], We have varied the heavy gauge boson mass between M G ut/2 < M x < 2M GGT , 
which reflects our ignorance of the GUT-scale mass spectrum. From this figure, we see 
that the existence of DM and Higgs multiplets produces a large effect on the proton decay 
lifetime. In particular, in the case of SU(4)cG»SU(2)iG)SU(2)^(8)Zl, the predicted lifetime 
is so small that the present bound has already excluded the possibility. This conclusion 
can be evaded, however, once the DM and i ? 2 Higgs multiplets are included in the theory 
as they raise the value of Mgut- Moreover, Model I is now being constrained by the proton 
decay experiments. In this case, the inclusion of the DM and Higgs multiplets decreases 
Mgut- Future proton decay experiments, such as the Hyper-Kamiokande experiment [42], 
may offer much improved sensitivities (by about an order of magnitude), with which we 
can probe a wide range of parameter space in both models. 

5.3 Non-equilibrium thermal dark matter 

Finally, we evaluate the relic abundance of DM produced by the NETDM mechanism [4] 
in Models I and II. In both of these models, the DM -0° is produced in the early Universe 
via the exchange of the intermediate-scale particles. Therefore, the production rate is 
extremely small and their self-annihilation can be neglected. In addition, the produced 
DM cannot be in the thermal bath since they have no renormalizable interactions with 
the SM particles. These two features characterize the NETDM mechanism; the DM is 
produced by SM particles in the thermal bath via the intermediate boson exchange, while 

17 For recent analyses on proton decay in SUSY GUTs, see Ref. [39]. 


21 




Figure 5: Proton lifetimes as functions of M X /M GUT . Blue solid and red solid lines 
represent Model I and Model II, respectively. Blue dashed and red dashed lines represent 
the cases for G int = SU(4 )c®SU(2) l ®SU(2) r andG mt = SU(A) C ® SU(2) L ®SU(2) R ® D 
when the DM and extra Higgs multiplets are not included. The shaded area shows the 
region which is excluded by the current experimental bound, r(p —> e + 7r°) > 1.4 x 10 34 years 

[40,41]. 


they do not annihilate with each other nor attain thermal equilibrium. I 11 what follows, we 
estimate the density of the DM produced via this mechanism and determine the reheating 
temperature which realizes the observed DM density. 

The Boltzmann equation for the DM if 0 is given by 


gD'dm 

dx 



Mpi 


Hy2 
x 2 eq ’ 


(28) 


with Ydm = it-bm/s and D e q = n eq /s, where «dm is the DM number density, n eq is the 
equilibrium number density of each individual initial state SM particle, and s is the 
entropy of the Universe; x = Mdm/F, with T being the temperature of the Universe; 
g* s and are the effective degrees of freedom for the entropy and energy density in the 

thermal bath, respectively; Mp\ = 1 /y/Grf — 1-22 x 10 19 GeV is the Planck mass; (av) is 
the thermally averaged total annihilation cross section of the initial SM particles, /, into 
the DM pair. When we derive Eq. (28), we neglect the DM self-annihilation contribution 
as discussed above. From now on, we assume g* s = g* p = g* for brevity. 

We evaluate the thermal averaged cross section (av) multiplied by the equilibrium 
number density squared n 2 q as 


(av)n 2 eq ~ 


T 


512vr 5 


'AM 2 

^ m dm 


dsfMiM^K^VS/T) Y. \M? , 


(29) 
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Figure 6: Diagram responsible for the DM production in Model II. 


where y/i denotes the center-of-mass energy, and K n {x) is the modified Bessel function 
of the second kind. Here, we have used the approximation mf •C V$i with rnj being 
the masses of the SM particles, since the particle production predominantly occurs at 
high temperature, and we have neglected the angular dependence of A4 for simplicity. In 
addition, we have assumed the initial particles follow a Maxwell-Boltzmann distribution 
and ignored statistical mechanical factors which may result from the Fermi-Dirac or Bose- 
Einstein distribution. |A4 | 1 indicates the sum of the squared amplitude over all possible 
incoming SM particles, as well as the spin of the final state. 

Next, we evaluate the amplitude A4 in each model. First, we consider the case of 
Model II. In this case, the dominant contribution comes from the tree-level Higgs-boson 
annihilation process displayed in Fig. 6. Here, if 0 , h , and 0 ° denote the DM, the SM 
Higgs boson, and the singlet component of the ( 15 , 1 , 1 )#, respectively, and the VEV ( 0 ) 
is given in Eq. (21). From the dimensional analysis, we estimate the contribution as 


\M\ 2 ~ c 


4-^dm 


ML 


(30) 


where c is a numerical factor which includes the unknown couplings appearing in the 
diagram. By substituting Eqs. (29) and (30) into Eq. (28), we have 


dY DM _ c /45\^ M p1 M dm 1 f 2 
dx ~ 1024tt 7 \ngj M 2 nt x 2 J 2x 


4a; 2 ) 2 Ki(t)dt . 


(31) 


When Mdm -C 7rh with Trh being the reheating temperature, the above equation is 
easily integrated to give 


(o) c f 45 \ 2 MpiT RH 
DM _ 64 it 7 (irg.) M' J a 


(32) 


where the superscript “(0)” implies the present-day value. On the other hand, the current 
value of is given by 


vHO) = ^ DM Per it 

DM m dmS (°) 


(33) 
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Figure 7: Examples of diagrams responsible for the DM production in Model I. 


where Odm is the DM density parameter and p^it is the critical density of the Universe. 
In the following calculation, we use Ddm ^ 2 = 0.12, p^J lt = 1.05 x 10 ~ 5 h 2 GeV • cm -3 , and 
s ( - 0 - ) = 2.89 x 10 3 cm -3 , with h the Hubble parameter. As a result, we obtain 


T rh ~ 2.7 x 10 4 GeV x 


01dm h 2 


3 

2 — 1 
9 *c 

10 4 


Mi 


DM 


100 GeV 


-l 


(34) 

This 


l3 -66 Q e y f rom the result in Table 6 . 

3 

is an unknown factor and 


0.12 

where we have set the value of M int = 10 

approximate formula is valid when Mdm •C Trh- Here, glc 
thus causes an uncertainty in the computation. For instance, if g* = 0(100) and the 

3 

quartic coupling of h and f> is ~ 0.3, then g*c~ l = 0(1O 4 ). Note that the perturbativity 
of the quartic coupling ensures that this factor cannot become too small. On the other 
hand, it also has an upper bound; if c is extremely small, then the one-loop gauge- 
boson exchange contribution dominates over the tree level. Taking this consideration into 

3 

account, we vary the value of g^M 1 by a factor of ten to estimate the uncertainty in the 
analysis given below. 

Next, we evaluate the relic abundance of the DM in Model I. In this case, there is 
no tree-level process for the DM production, since the DM does not couple to the singlet 
component </>° in the (1,1, 3 )r. Therefore, the DM is produced at the loop level. In Fig. 7, 
we show examples of one-loop diagrams which give the dominant contribution to the DM 
production. The amplitude is then estimated as 


\ I Ad 1 2 — ——— 

1 (16vr 2 ) 2 


4^dm 


Ml 


(35) 


int 


where we have included the one-loop factor. After a similar computation, we obtain 


yd 0 ) ~ 
Dm - 


/ 45 \ 2 M Pl T } 


- 


647 t 7 (167 t 2 ) 2 \ng 


RH 


ML 


and 


Trh — 4.6 GeV x 


f^DM ^ 2 

0.12 


g*c' 1 
10 5 


M 


DM 


GeV 


-i 


(36) 


(37) 
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(a) Model I 


(b) Model II 


Figure 8: Reheating temperature as a function of DM mass. Pink band shows the theoret¬ 
ical uncertainty. 


on the assumption of Mdm •C T rh . Here, we have set M int = 10 8 ' 08 GeV. 

In Fig. 8, we plot the predicted reheating temperature as a function of the DM mass 
after numerically integrating Eq. (31). The left and right panels show the cases of Model 
I and II, respectively. The pink band shows the uncertainty of the calculation, which 
we estimate by varying the unknown factor by a factor of ten. It turns out that when 
Mdm Trh, hr the case of Model I, only a small DM mass is allowed and the reheating 
temperature must be quite low. In Model II, on the other hand, DM with a mass of 
around the electroweak scale accounts for the observed DM density with an acceptably 
high reheating temperature. For a larger Mdm, hr both models, the DM relic abundance 
can only be explained in the narrow strip region where Mdm — Trh. 

6 Lonely Singlet Fermion Dark Matter 

In the above discussion, we have assumed that there exists a DM multiplet (as well as 
extra Higgs multiplets) above the intermediate scale, and studied how the presence of the 
additional fields affect the gauge coupling running in such models. As seen in Sec. 3.2, 
these fields can indeed improve the solutions for both the intermediate and GUT scales, 
which allow the models to evade the limit from the proton decay experiment and to explain 
light neutrino masses via the seesaw mechanism. Before concluding our discussion, we 
briefly consider another possibility in this section; that is, we have only a singlet DM 
fermion on top of the standard SO(IO) setup discussed in Sec. 3.1. In this case, the 
DM, of course, cannot affect the gauge coupling running, and thus it does not solve 
the problems regarding the low intermediate/GUT scales in the ordinary SO(IO) GUT 
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models. Since there may be another solution to these problems, it is worthwhile studying 
this possibility as well. 

In fact, we can easily construct such a model by exploiting an appropriate Higgs 
held at the GUT scale and fine-tuning its VEV so that only the singlet fermion DM 
has a mass much lighter than the GUT scale. For example, let us consider the case of 
SU(4)c <8> SU(2)i <8> SU(2) fi <g) D. In this case, the singlet held under the intermediate 
gauge interactions, ( 1 , 1 , 1 ), is contained in a 54 or 210 of SO(IO). Since only the 210 
can have a Yukawa coupling to the R± = 54^ Higgs held, we focus on the case where the 
singlet DM fermion is a component of the 210 held. In this case, both Majorana and 
Dirac fermions can couple to the R\ Higgs. Then, by hne-tuning the Yukawa coupling, we 
can make only the singlet component have a light mass, as is done in Sec. 4. Similarly, we 
can obtain other models with different intermediate gauge groups by using appropriate 
multiplets for the helds which contain the singlet DM. 

The NETDM mechanism again works for this singlet DM through the R\ Higgs ex¬ 
change process at tree level, with a diagram similar to that illustrated in Fig. 6. Following 
the discussion given in Sec. 5.3, we can readily evaluate the reheating temperature required 
to produce the right amount of DM. When Mdm Trh, we have 


Trh — 1.3 x 10® GeV x 


Odm h 2 \fg*c 1N \ / Mdm \ ' / Mgut "V 
0.12 ) V 10 4 ) V100 GeV ) V10 16 GeV ) 


(38) 


Compared with Model I and II, the present scenario in general predicts a high reheating 
temperature, as the production occurs via the GUT-scale particle exchange. Such a high 
reheating temperature may be consistent with thermal leptogenesis [43]. 

As for proton decay and neutrino masses, the consequence of the singlet DM models 
is the same as that without DM. Thus, for G- mt = SU(4)c 8> SU(2)^ (8) SU(2) R <g> D and 
SU(4) C (8> SU(2 ) l ® U(1) r , the proton decay constraints are still problematic, and thus it 
may be required that we assume a relatively heavy GUT-scale gauge boson when compared 
to the GUT scale. Other intermediate groups are not suitable for the explanation of 
neutrino masses. The solution discussed in Sec. 5.1 can again be exploited in these cases. 


7 Conclusion and discussion 

For over 40 years now, we have wondered whether grand unification is actually realized 
in nature. Its simplicity, its capacity for an explanation of charge quantization and the 
apparent focusing of the gauge couplings as they run to high energy has kept grand 
unification (supersymmetric or not) at the center of most ultra-violet completions of the 
SM, though experimental verification is still lacking. 

On the other hand, we know from the existence of neutrino masses, the baryon asym¬ 
metry of the Universe and the existence of DM that there must be new physics beyond 
the SM. The presence of a natural DM candidate in SUSY extensions of the SM (with 
conserved .R-parity) is often taken to be one of the motivations for low-energy SUSY. 
The ingredients for the baryon asymmetry are contained in most grand unified theories 
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(supersymmetric or not) including SU(5) and SO(IO), and while a neutrino seesaw can 
be accomplished in SU(5) (by including the right-handed neutrino as a SU(5) singlet), it 
is more natural in SO(IO). 

We have, here, examined several breaking schemes of SO(IO) which lead to gauge 
coupling unification (by altering the SM running of the gauge couplings at an intermediate 
scale), and contain a remnant symmetry which can account for the stability of DM. 
Having established the possible intermediate-scale gauge groups capable of both gauge 
coupling unification and of supporting a stable DM candidate, we considered specific 
possible representations (of dimension no larger than 210 for simplicity) which contain a 
suitable non-degenerate SM singlet DM candidate. If the DM candidate couples to the 
SM only through intermediate scale fields, it may never equilibrate in the early Universe 
after reheating, and its production from the thermal bath is an example of the NETDM 
scenario. Despite the fact that there are several possible intermediate-scale gauge groups 
to consider and many possible representations for the DM candidate and intermediate- 
scale Higgs fields needed to break the degeneracy in the DM multiplet, we found only two 
surviving models: one each based on SU(4)c <g) SU(2)£, g) SU(2) R and SU(4)c <S> SU(2 ) L g 
SU(2)i?<g)D, with DM contained in a (1,1, 3)b G 45b and (15,1,1 )w G 45 w, respectively. 

Both of the surviving models are capable of producing light neutrino masses (though 
it is more difficult in Model I due to its relatively low intermediate scale). We also showed 
that while the proton decay lifetime (to e + 7T°) is at least a factor of two longer than 
the current experimental bound for Mx/Mqut > 1/2 in Model I, the current bound 
excludes masses M x /Mq VT < 0.7, and higher masses may be probed in future proton 
decay experiments. Finally, within the NETDM production scenario, we have related our 
two models to a specific reheat temperature after inflation needed to obtain the current 
relic density. While Model II predicts a reheat temperature which easily allows for (non 
thermal) leptogenesis [43,44], the reheat temperature in Model I is rather low and presents 
a challenge for baryogenesis. 
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Appendix 


A Input parameters 


The values for the input parameters we have used in this paper are summarized in Table 7. 
They are taken from Ref. [45] except for the top-quark pole mass and the Higgs mass, for 
which we use the values given in Refs. [46] and [47], respectively. In this table, the gauge 
coupling constants are defined in the MS scheme, and thus we convert them to the DR 
scheme at the electroweak scale using the one-loop relation [48]: 

<7a( m .z) D R — gaymzjMS ( 1 + 24^ ) > ( 39 ) 

where C(G a ) is the quadratic Casimir invariant. For the mass of the top quark, we convert 
the pole mass to its MS mass by using [45] 

(40) 


m 


MS 


(m^ s ) = m t (1 — 


4 a s (m\ 
37r 


from which we obtain the MS top Yukawa coupling, 
given by 


% DR = 


i + 


OL i 


+ 


3ct2 


4807T 327T 


The DR Yukawa 


Q?3 

37T 


coupling is then 


(41) 


Table 7: Input parameters [45-47], 


Strong coupling constant 

a s (m z ) 

0.1185(6) 

QED coupling constant 

a{m z ) 

1/127.944(14) 

Fermi coupling constant 

G f 

1.1663787(6) x 10~ 5 GeV“ 2 

Weak-mixing angle 

sin 2 9 w (m z ) 

0.23126(5) 

Z-boson mass 

m z 

91.1876(21) GeV 

Top pole mass 

m t 

173.34(82) GeV 

Higgs mass 

m h 

125.15(24) GeV 


B Renormalization group equations 


In this section, we summarize the RGEs and the matching conditions used in text. The 
two-loop RGEs [49] of the gauge coupling constants g a are written as 

1.(1) 3 r 3 


dg a _ 

^ dpL 167r 2i, “ 


9a + 


9 a 


(167T : 


2)2 


E b 

■ 6=1 


( 2 ) 2 
ab 9b 


CalJt 


(42) 



Below, we will give the coefficients in each theory discussed in this paper. For the con¬ 
tribution of Yukawa couplings, we include them only in the SM running, as unknown 
Yukawa couplings appear above the intermediate scale. Their effects should be taken into 
account as theoretical uncertainties. All of the one-loop RGEs have been checked with 
the code PyR@TE [50], and more importantly, the two-loop RGEs have been computed 
with this code. 


B.l Standard Model 

In the SM, we have 



( 41/10 ^ 
-19/6 

b {2) - 

% U ab — 

A99/50 27/10 44/5^ 
9/10 35/6 12 

1 Ca 

( 17/10^ 

3/2 


l " 7 ) 


^11/10 9/2 -26 j 


l 2 ) 


Here, a = 1,2,3 correspond to U(l), SU(2) L , and SU(3)c, respectively, with the U(l) 
gauge coupling constant normalized as gi = \j3/3g'. Since the top Yukawa coupling 
contributes to the running of the gauge couplings at the two-loop level, it is sufficient 
to consider the one-loop RGE for the top Yukawa coupling. Furthermore, we can safely 
neglect the contribution of the other Yukawa couplings. Thus, the relevant RGE is 


g-rvt = 

an I07 t z 


9 2 17 2 9 2 2 

2 Vt ~ 2 q . 9 i ~ 4^2 - 


(44) 


B.2 SU(4) C 0 SU(2) l 0 SU(2 ) R 

As discussed in Sec. 3.1, above the intermediate mass scale, the theory contains the SM 
fermions, the gauge bosons, the (10, l,3)c held, and the (l,2,2)c Higgs field. The 
beta-function coefficients in this case are given by 


b m = 

( _3 l 

11/3 

b (2) - 

J U ab ~ 

( 8 3 45/2 ^ 

3 584/3 765/2 


{- 23 / 3 ) 


^9/2 153/2 643/6y 


(45) 


where a = 2 L, 2R, 4 correspond to SU(2 )l, SU(2)#, and SU(4)c, respectively. The match¬ 
ing conditions at the intermediate mass scale are 

1 _ 3 1 2 1 

5 gl R {M- mt ) + 5 gj(M int ) 

<72 ( A/]i n t) = 9 2L(M int ) , 

int) = g 4 (Mi nt ) . (46) 
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B.3 SU(4)c 0 SU(2 ) l <8) SU(2) i? 0 D 

In this case, the (10, 3,1 )c field is added to the previous theory. The beta-function 
coefficients then become 

( 11/3 \ /584/3 3 765/2 ^ 

b { a ] = 11/3 , b$ = 3 584/3 765/2 , (47) 

^—14/3 y yl53/2 153/2 1759/6, 

where a = 2L,2i?, 4 correspond to SU(2) L , SU(2)/j, and SU(4)<®, respectively. 

B.4 SU(4) c ®SU(2) L <g>U(l) ii 

This theory contains the SM fermions, the gauge bosons, the (10,1, l)c field, and the 
(1, 2, |) Higgs field. The beta-function coefficients in this case are given by 

/—19/6\ /35/6 1/2 45/2 N 

b ( a ] = 15/2 , b$= 3/2 87/2 405/2 , (48) 

v -29/3 ) y 9/2 27/2 -101/6, 

where a = 2 L, 1R, 4 correspond to SU(2 )l, U(1)r, and SU(4)c, respectively. The match¬ 
ing conditions at the intermediate mass scale are 

1 _ 3 1 2 1 
gf(M int ) ~ 5 gf R (M int ) + 5 gj(M int ) ’ 

<72(44i nt ) = g 2L (M int ) , 

fl's(44i n t) = . (49) 

B.5 SU(3)c 0 SU(2 ) l 0 SU(2) i? 0 U(1 ) B -l 

This theory contains the SM fermions, the gauge bosons, the (1,1, 3, —2 )c field, and the 
(1, 2, 2,0) Higgs field. The beta-function coefficients in this case are given by 

/ —3 \ / 8 3 3/2 12 

6 (i )= -V 3 h {2 )= 3 80/3 27/2 12 

“ 11/2 ’ ab 9/2 81/2 61/2 4 

\ —7 / \9/2 9/2 1/2 -26 

where a = 2L,2R, BL,3 correspond to SU(2) L , SU(2) R , U(1 )b_l and SU(3)c, respec¬ 
tively. The U(1 )b-l charge is normalized such that it satisfies the normalization condi¬ 
tion of the SO(10) generators: Tb~l = >/3/8 [B — L). The matching conditions at the 
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intermediate mass scale are 


1 _ 3 1 2 1 
0 i(-Mmt) 5 <4(M int ) + 5 g 2 BL (M int ) ’ 

02(M nt ) = g2L(M int ) , 

03 (-Mint) = 03(Mnt) • (51) 


B.6 SU(3) C 0 SU(2) l 0 SU(2) i? 0 U(1) B _ L 0 D 

For this left-right symmetric theory, the (1,3,1,2)^ held is added to the previous case. 
The beta-function coefficients are then modified to 


where a 
tivcly. 
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SU(2) r 

, U(1) B _ 
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(52) 


B.7 SU(3) c 0SU(2) t ®U(l) R ®U(l) iJ _ £ 

This theory contains the SM fermions, the gauge bosons, the (1,1,1, —2 )c field, and the 
(1, 2, 1/2, 0) Higgs field. The beta-function coefficients in this case are given by 


/-19/6\ 


/ 35/6 
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15/2 
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3/2 

1/2 

—26/ 


(53) 


where a = 2 L, 1 R, BL , 3 correspond to SU( 2 ) L , U(l)#, U( 1 )b_^ and SU(3)c, respectively. 
The matching conditions at the intermediate mass scale are 


1 _ 3 1 2 1 
0 i 2 (M int ) 5 g 2 1R (M int ) + 5 g 2 BL (M iDt ) ’ 

0 2 (M int ) = g2L(M int ) , 

03 (-Mint) = 03 (-Mint) • (54) 


B.8 Model I 

For DM model I, a (1,1, 3 )d Dirac fermion and a (1,1, 3)r real scalar field are added to 
the theory described in Appendix B.2. The beta-function coefficients are then computed 
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( -3 N 
20/3 

b (2) - 

1 U ab ~ 

( 8 3 45/2 ^ 

3 740/3 765/2 


y—23/3y 


y9/2 153/2 643/6y 


where a = 2L,2R, 4 correspond to SU(2 )l, SU(2)r, and SU(4) G , respectively. 


(55) 


B.9 Model II 


For DM model II, a (15,1,1 )w Weyl fermion and a (15,1, 1 )r real scalar field are added to 
the theory described in Appendix B.3. The beta-function coefficients are then computed 
as 



( H/3 ^ 
11/3 , 


v - 4 / 3 ! 


/584/3 
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765/2 ^ 
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584/3 

765/2 

, (56) 

v 153/2 
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where a = 2L,2R, 4 correspond to SU(2) L , SU(2)^, and SU(4) G , respectively. 


C One-loop formulae for gauge coupling unification 


At the one-loop level, the gauge coupling RGEs are easily solved analytically. By using 
the solutions, we can obtain analytic expressions for M int , M GUT , and «gut as follows: 


Mint = rn z exp 


M gut = rriz exp 


27r(h x n ) • «_! 

(b x n) ■ b 

27t(A b x n) • «_i 

(b x n) • b 
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(60) 


where A b = b — b, and b a and b a denote the beta-function coefficients below and above 
the intermediate scale, respectively. The U(l) beta function above the intermediate scale 
is given by a linear combination of the beta functions of the intermediate gauge group. 
For instance, in the case of SU(4) G <g) SU(2)i g) SU(2) fi , we have 


I 2 , 3 , 

b\ — -»4 + 

5 5 


(61) 
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Similar expressions are obtained for other intermediate groups. Notice that the compo¬ 
nents of the beta-function coefficients which are proportional to n do not affect Mqut 
and Mi nt , as one can see from the formulae. Therefore, if one adds a multiplet to, e.g., the 
SU(4)c <E> SU(2 ) l (g) SU(2 ) r theory whose contribution to the beta-function coefficients is 
A5 4 = A b 2L = A b 2R , then the multiplet does not alter M GUT and M int at the one-loop 
level. 

We also note that physics above the intermediate scale gives negligible effects on the 
determination of M- mt in the presence of the left-right symmetry. We can see this feature 
by using Eq. (57). Let us consider the case of SU(4 )c<8>SU(2)l(8)SU(2)b<8) D. In the left- 
right symmetric theories, the beta functions of the SU(2) L and SU(2) R gauge couplings 
should be the same. Therefore, we have b 2L = b 2R , and 

b x n = ( b 2L - b 4 )c , (62) 

with 

(63) 

Therefore, Eq. (57) reads 

T 27 TC • ct_ il . 

Mint = mz exp --— , (64) 

c • b 

and thus, the intermediate scale does not depend on the beta function above M int . One 
can also see this feature by noting that above the intermediate scale g 2 L = g 2 n holds 
at any scale. Hence, the intermediate scale corresponds to a point at which g 2 l becomes 
equivalent to g 2 R, which is determined only by the running below Mi nt . A similar argument 
holds in the case of SU(3)c ® SU(2 )l ® SU(2)^ ® U(1)b_l ® D. 

D Proton decay in SO(IO) —> SU(4) 0 SU(2) 0 SU(2) 

Here, we give details of the calculation for the proton decay lifetime in the intermediate- 
scale scenario. We consider the case of SO(IO) —» SU(4) ® SU(2) ®) SU(2), which was 
discussed in Sec. 5.2. 

In non-SUSY GUTs, proton decay is induced by gauge interactions. The relevant 
interactions are written as 

Cm, = ^ [(Q),.,.X“'-Pr(£ c ). + (Q)a,f“ r P l(L c ) t + + h.c.] , 

(65) 

where 
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and X denotes the superheavy gauge bosons which induce the baryon-number violating 
interactions; (?gut is the unified gauge coupling constant; a,b,c, i,j, and r,s are the 
SU(3)c, SU(2 )l, and SU(2)# indices, respectively; Pr/l = (1 ±75)/2 are the chirality 
projection operators. 

After integrating out the SO(IO) gauge fields X, we obtain the dimension-six proton 
decay operator. The operator is expressed in a form that respects the intermediate gauge 
symmetry, SU(4) SU(2) ® SU(2): 

£ eff = C(M G xjt) ■ e i3 e rs e a ^ 5 {W) ai P L ^(Wy r P R ^ s , (67) 

where a, (3,... denote the SU(4) indices, and T is given in Eq. (13). Notice that 

e ij e kl e a frs(* 5 ) ai PL* f>i (&)' yk PL* 51 = e rs e tu e a ^ 5 {W) ar P R ^ S {Wy P R ¥ U = 0 , ( 68 ) 


and thus the operator in Eq. (67) is the unique choice. At tree level, the coefficient of the 
effective operator is evaluated as 


C(M gut) 


ffGUT 

2 M\ ’ 


(69) 


with Mx the mass of the heavy gauge field X. Here, we have neglected fermion flavor 
mixings [51] for simplicity. 

The Wilson coefficient is evolved down to the intermediate scale using the RGE. The 
renormalization factor is computed to be [52] 




0(4 ( A(/int) 

4 b 4 

0)2 L(Mint) 

4 h l 

0)2 J?(Mi nt ) 

ckgut 


CtGUT 


ckgut 


4^2 R 

C(^gut) • 


(70) 


At the intermediate scale, the SU(4) <® SU(2) ®> SU(2) theory is matched onto the SM. 
The effective Lagrangian is written as 


C eS = J2 C I°Ir (71) 

1=1 

with the effective operators given by [53-55] 

o, = 4) , 

0 2 = €^ £y (QrQ?)(« c R e K ) , 

o 3 = e^,MQTQi k )(QL L i) , 

— e a bc(Uftd b R )(u c R e R ) ■ (72) 

We evaluate the coefficients Cj as 

Cl (Mint) = 4 C(M int ) , 

C 2 (M int ) = —4(7(Mi nt ) , 

C 3 (M int ) = C 4 (M int ) = 0. (73) 
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We then run down the coefficients to the clectroweak scale. The renormalization 
factors are given by [55] 


Ci(mz) 

C 2 (m z ) 


' a 3 (m z ) " 

_2_ 

b 3 

’ oi 2 {m z ) ' 

9 

4b 2 

' ai (m z ) ' 

11 

20 b 1 

Cl (Mint 

_a 3 (M int )_ 


_a 2 (Adi nt ) _ 


_ai(M int ) _ 

a 3 (m z ) 

2 

b 3 

at 2 (m z ) 

9 

4b 2 

ot\{m z ) 

23 

20b! 

C2(M int 

a 3 (M int )_ 


_a 2 (Mi nt ) _ 


_o;i(M'i n t) _ 


(74) 

(75) 


Note that the beta-function coefficients should be appropriately modified when the num¬ 
ber of quark flavors changes. Below the electroweak scale, the QCD corrections are the 
dominant contribution. By using the two-loop RGE given in Ref. [56], we compute the 
Wilson coefficients at the hadronic scale ythad as 


Ci{n 


had / 


(/^had) 

6 

25 

a s (m b ) 

6 

23 

®s(R had) + 

173 

825 

~a s (m b ) + ^' 

_ a s (m b ) _ 


_a s (m z )_ 


_ a s (m b ) + - 


a s (m z ) + p|Q 


430 

2001 


Ci{m z ) 




(76) 


with i = 1,2. 

In non-SUSY GUTs, the dominant decay mode of proton is p —> 7r°e + . The partial 
decay width of the mode is computed as 


T(p —> 7r°e + ) 



(77) 


where rn p and m*- are the masses of the proton and the neutral pion, respectively, and 


A L = Ci(n had )(iT 0 \(ud) R u L \p) , 

Ar = 2C 2 (^had) (7T° | (ltd) lUr \p) . (78) 


The hadron matrix elements are evaluated with the lattice QCD simulations in Ref. [57]. 
We have 


(n°\(ud) R u L \p) = (n 0 \(ud) L u R \p ) = -0.103(23)(34) GeV 2 , (79) 

with /ih a d = 2 GeV. Here, the hrst and second parentheses indicate statistical and sys¬ 
tematic errors, respectively. 
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